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Abstract 

We present a translationally invariant formulation of the no-core shell model 
approach for few-nucleon systems. We discuss a general method of antisym- 
metrization of the harmonic-oscillator (HO) basis depending on Jacobi co- 
ordinates. The use of a translationally invariant basis allows us to employ 
larger model spaces than in traditional shell-model calculations. Moreover, 
in addition to two-body effective interactions, three- or higher-body effective 
interactions as well as real three-body interactions can be utilized. In the 
present study we apply the formalism to solve three and four nucleon sys- 
tems interacting by the CD-Bonn nucleon-nucleon (NN) potential in model 
spaces that include up to 34:hQ and 16h~Q HO excitations, respectively. Re- 
sults of ground-state as well as excited-state energies, rms radii and magnetic 
moments are discussed. In addition, we compare charge form factor results 
obtained using the CD-Bonn and Argonne V8' NN potentials. 



I. INTRODUCTION 

Various methods have been used to solve the few-nucleon problem in the past. The 
Faddeev method |j| has been successfully applied to solve the three-nucleon bound-state 
problem for different nucleon-nucleon potentials |@-fll- For the solution of the four- nucleon 
problem one can employ Yakubovsky's generalization of the Faddeev formalism || as done, 
e.g., in Refs. f| or |7]]. Alternatively, other methods have also been succesfuliy used in the 
past, such as, the correlated hyperspherical harmonics expansion method j8],|9| or the Green's 
function Monte Carlo method ||10|| . Apart from the coupled cluster method fTT| , |T2|j appli- 



cable typically to closed-shell nuclei, the Green's function Monte Carlo method is the only 
approach at the present time, for which exact solutions of systems with A > 4 interacting 
by realistic potentials can be obtained, with the current limit being A = 8. 
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On the other hand, when studying the properties of more complex nuclei, one typically 
resorts to the shell model. In that approach, the single-particle harmonic-oscillator basis 
is used and the calculations are performed in a truncated model space. Instead of the free 
NN potential, one utilizes effective interactions appropriate for the truncated model space. 
Examples of such calculations are the large-basis no-core shell-model calculations that have 
recently been performed [|I^,|14[] . In these calculations, the effective interaction is determined 
for a system of two nucleons bound in a HO well and interacting by the nucleon-nucleon 
potential. We note that the use of a HO basis is crucial for insuring that the center-of-mass 
motion of the nucleus does not mix with the internal motion of the nucleons. On the one 
hand, this approach is limited by the model-space size, and, on the other hand, by the fact 
that only a two-body effective interaction is used despite the fact that higher-body effective 
interactions might not be negligible. 

It is possible, however, to re-formulate the shell-model problem in a translationally- 
invariant way Recently, we combined the no-core shell- model approach to the three- 

and four-nucleon systems with the use of antisymmetrized translationally invariant HO basis 



18|,p79[| . That allowed us, due to the omission of the center of mass, to extend the shell-model 



calculations to model spaces of 32hQ and 14hQ excitations above the unperturbed ground 
state for the A = 3 and A = 4 systems, respectively. In addition, that approach made it 
possible to employ the three-body effective interactions in the A = 4 calculations. 

In the present paper we simplify and generalize this approach so that it is applicable 
to an arbitrary number of nucleons. In particular, we discuss in detail how to construct 
an antisymmetrized HO basis depending on Jacobi coordinates. We present an iterative 
formula for computing the antisymmetrized basis for A nucleons from the antisymmetrized 
basis for A — 1 nucleons. Further, we discuss how to transform the antisymmetrized states 
to bases containing different antisymmetrized subclusters of nucleons. 

We also describe the effective interaction derivation from a different perspective than in 
our previous papers. Namely, we point out the connections between the no-core shell-model 
approach and the unitary- model-operator approach |20| . Let us remark that a fundamental 
feature of the effective interactions that we employ is the fact that with the increasing 
model-space size the effective interactions approach the bare NN interaction. Therefore, in 
principle, our approach converges to the exact few-nucleon solution. 

The basic advantage of this formalism is, first, the fact that larger model spaces can be 
utilized than in the standard shell-model calculations, because the center-of-mass degrees of 
freedom are omitted and because a coupled basis with good J and T is used. Second, due 
to the flexibility of HO states depending on Jacobi coordinates, different recouplings of the 
basis are possible. Consequently, not only can two-body effective interactions be utilized, 
but also three- or higher-body effective interactions as well as real three-body interactions. 
On the other hand, because the antisymmetrization procedure is computationally involved, 
the practical applicability of the formalism is limited to light nuclei. In the present formula- 
tion we expect that significant improvement over the traditional shell-model results can be 
achieved for A < 6. 

We apply the formalism to solve three- and four-nucleon systems interacting by the non- 
local momentum-space CD-Bonn NN potential [^TJ. The present calculations are done in 
larger model spaces than those used in Refs. |18|JT9|] . Also, there have not been any published 
results so far for the A = 4 system interacting by the CD-Bonn NN potential. In addition 
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to the calculation of ground-state and excited-state energies, point-nucleon rms radii and 
magnetic moments, we evaluate electromagnetic (EM) and strangeness form factors in the 
impulse approximation. We plan to present the results of calculations for A = 5 and A = 6, 
using the present formalism, separately. 

In Section O, we discuss the standard no-core shell-model formulation, i.e., the Hamilto- 
nian and effective interaction calculation. The construction of the translationally invariant 
HO basis is described in Section |mi Results for A = 3 and A = 4 systems interacting by 
the CD-Bonn NN potential are given in Section |IV|. In Section [V|, we present concluding 
remarks. 



II. NO-CORE SHELL-MODEL APPROACH 



A. Hamiltonian 



In the no-core shell-model approach we start from the one- plus two-body Hamiltonian 
for the yl-nucleon system, i.e., 



E 

i<j=i 



(i) 



where m is the nucleon mass and Vt 



N{ri 



the NN interaction. In order to simplify the 
notation, the spin and isospin dependence is omitted in the interaction term in Eq. (JI|). We 



can use both coordinate-space dependent NN potentials, such as the Reid, Nijmegen [22] or 



Argonne |TDf as well as momentum-space dependent NN potentials, such as the CD-Bonn 
|21|| . In the next step we modify the Hamiltonian (||) by adding to it the center-of-mass HO 
potential ^Am£l 2 R 2 , R = i^fj. This potential does not influence intrinsic properties 
of the many-body system. It allows us, however, to work with a convenient HO basis and 
provides a mean field that facilitates the calculation of effective interactions. The modified 
Hamiltonian, depending on the HO frequency Q, can be cast into the form 



E 

i=l 



Pi i 1 r>2 -Q. 

— m\ Ir? 

2m 2 



A 

+ E 

i<j=l L 



rriQ 2 
~2A 



(2) 



The interaction term of the Hamiltonian (H) depends only on the relative coordinates. The 
one-body term in Eq. (|2|) can be re-written as a sum of the center-of-mass term, and a term 
depending on the relative coordinates. 

The shell-model calculations are performed in a finite model space. Therefore, the in- 
teraction term in Eq. (fj) must be replaced by an effective interaction. In general, for an 
A-nucleon system, an A-body effective interaction is needed. In practice, the effective in- 
teraction is usually approximated by a two-body effective interaction. In the present study 
we will also employ a three-body effective interaction. As approximations are involved in 
the effective interaction treatment, large model spaces are desirable. In that case, the cal- 
culation should be less affected by any imprecision of the effective interaction. The same 
is true for the evaluation of any observable characterized by an operator. In the model 
space, renormalized effective operators are required. The larger the model space, the less 
renormalization is needed. 
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As the Hamiltonian (@) differs from the Hamiltonian Ha ([]]) only by a center-of-mass 
dependent term, no dependence on Q should exist for the intrinsic properties of the nucleus. 
However, because of the approximations involved in the effective interaction derivation, a 
dependence on Q appears in our calculations. This dependence decreases as the size of the 
model-space is increased. 



B. Effective interaction theory in Lee-Suzuki approach 



In our approach we employ the Lee-Suzuki similarity transformation method p3| , |24 
which yields a starting-energy independent Hermitian effective interaction. In this subsection 
we recapitulate general formulation and basic results of this method. Applications of this 
method for computation of two- or three-body effective interactions are described in the 
following subsections. 

Let us consider an arbitrary Hamiltonian H with the eigensystem E^, \ k), i.e., 

H\k) = E k \k) . (3) 

Let us further divide the full space into the model space defined by a projector P and the 
complementary space defined by a projector Q, P + Q = 1. A similarity transformation of 
the Hamiltonian e'^He^ can be introduced with a transformation operator u satisfying the 
condition uo = QujP. The transformation operator is then determined from the requirement 
of decoupling of the Q-space and the model space as follows 

Qe^He^P = . (4) 

If we denote the model space basis states as \ap), and those which belong to the Q-space, 
as \cxq), then the relation Qe~ u He u P\k) = 0, following from Eq. (01), will be satisfied for a 
particular eigenvector \k) of the Hamiltonian (|3|), if its Q-space components can be expressed 
as a combination of its P-space components with the help of the transformation operator ou, 
i.e., 

(a Q \k) = Y,(aQ\uj\a P }(ap\k) . (5) 

If the dimension of the model space is dp, we may choose a set /C of dp eigenevectors, for 
which the relation (|5]) will be satisfied. Under the condition that the dp x dp matrix (otp\k) 
for \k) G K, is invertible, the operator u can be determined from (|^) as 

(a Q \u\a P ) = J2( a Q\ k )(k\ a p) > ( 6 ) 

fcG/C 

where we denote by tilde the inverted matrix of (ap\k), e.g., J2 ap {k\ a p){ a pW) = &k,k'i f° r 
k, k! eJC. 

The Hermitian effective Hamiltonian defined on the model space P is then given by |24 



H, 



cff 



P{l+Ju)P 1/2 PH{P + QuP) P{l + Ju)P (7) 



-1/2 
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By making use of the properties of the operator u, the effective Hamiltonian H e g can be 
rewritten in an explicitly Hermitian form as 



H, 



eft 



-1/2 



P(l + Ju)P (P + Pu ] Q)H(QuP + P) 



P(l 



^lj)P 



-1/2 



(8) 



With the help of the solution for uj (|6]) we obtain a simple expression for the matrix elements 
of the effective Hamiltonian 

(a P \H eS \a P >) = E E E (^K 1 + ^^y 1/2 \(yp")(a pi '\k)E k (k\a p/// ) 

fcG/C «p» OLpin 

x(a P /»|(l + JLu)- 1/2 \a p/ ) . (9) 
For computation of the matrix elements of (1 + cu^cu)^ 1 ^ 2 , we can use the relation 

(a P \ (1 + Jw) \a P ») = E (a P \k) (k\a P „) , (10) 

fcG/C 

to remove the summation over the Q-space basis states. The effective Hamiltonian @ 
reproduces the eigenenergies Ek, k G K in the model space. We note that the relation (|9p 
used together with (pTOj ) does not contain any summation over the Q-space basis and, thus, 
represents a simplification compared to the formulas we presented in our previous papers 
TjUTBUBfl, though it is fully equivalent. 



C. Unitary transformation of the Hamiltonian and the two-body effective interaction 

We now return to our problem, namely, to derive the effective interaction corresponding 
to a chosen model space for the 'particular Hamiltonian (|2|). Let us write the Hamiltonian 
(0) schematically as 

A A 

flS = E^+E^- ( n ) 

i=l i<j=l 



According to Providencia and Shakin pSj , a unitary transformation of the Hamiltonian, 
which is able to accomodate the short-range two-body correlations in a nucleus, can be 
introduced by choosing a two-body antihermitian operator S, such that 

H = e~ s H%e s . (12) 

Consequently, the transformed Hamiltonian can be expanded in a cluster expansion 

n = n^ + n {2) + n {3) + ... , (13) 

where the one-body, two-body and three-body terms are given as 

H (1) = I>, ( 14a ) 

i=l 

= E % > (i4b) 

i<j=i 

A 

^ (3) = E ( i4c ) 

i<j<k=l 
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with 



V 12 = e-^ih + h 2 + V 12 )e Sw - (/i! + h 2 ) , (15a) 
V123 = e- Sl23 (/n + h 2 + h 3 + V 12 + V 13 + V 23 )e s ™ 

-(ht + h 2 + h 3 + V 12 + V 13 + V 23 ) , (15b) 

and 5i23 = 5*12 + 5*23 + S 3 \. In the above equations, it has been assumed that the basis states 
are eigenstates of the one-body, in our case HO, Hamiltonian J2t=i hi- 

If the full space is divided into a model space and a Q-space, using the projectors P and 
Q with P + Q = 1, it is possible to determine the transformation operator Si 2 from the 
decoupling condition 

Q 2 e- Sl *{h 1 + h 2 + V 12 )e Sl *P 2 = , (16) 

and the simultaneous restrictions P 2 S\ 2 P 2 = Q 2 S\ 2 Q 2 = 0. Note that two-nucleon-state 
projectors appear in Eq. (|T6"D, whose definitions follow from the definitions of the A-nucleon 
projectors P, Q. This approach, introduced by Suzuki and Okamoto and referred to as the 
unitary- model-operator approach (UMOA) PUj , has a solution that can be expressed in the 
following form 

5*12 = arctanh(a; — uy) , (17) 

with the operator u satisfying to = Q 2 uP 2 . This is the same operator that solves Eq. (f|) 
with H = h 1 + h 2 + Viz. Moreover, P 2 e" Sl2 (/ii + h 2 + V u )e Sl2 P 2 is given by Eq. ([|) with 
H = h 1 + h 2 + V 12 . 



D. Two-body effective interaction calculation 



We compute the two-body effective interaction according to Eq. (|15a|) using the decou- 
pling condition ([IB]). For the two-nucleon Hamiltonian we employ 

H? = H 02 + V 12 = ^ + + V N (V2f) - , (18) 

where r = \[\{f\ — r 2 ) and p = \J\(j>\ — P2) an d where Hq 2 differs from h\ + h 2 by the 
omission of the center-of-mass HO term of nucleons 1 and 2. Our calculations start with 
exact solutions of the Hamiltonian ([TJ) and, consequently, we construct the operator uj and, 
then, the effective interaction directly from these solutions by application of the relations 
(|) and (|) with E k , \ k) obtained from the solution of the Schrodinger equation (|3|) for 
the Hamiltonian H = H 2 . The relative-coordinate two-nucleon HO states used in the 
calculation are characterized by quantum numbers \nlsjt) with the radial and orbital HO 
quantum numbers corresponding to coordinate r and momentum p. Typically, we solve the 
two-nucleon Hamiltonian (p!8[) for all two-nucleon channels up to j = 6. For the channels 
with higher j only the kinetic-energy term is used in the many-nucleon calculation. The 
model space is defined by the maximal number of allowed HO excitations iV max from the 
condition 2n + I < iV max . 
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In order to construct the operator uj (J6|) we need to select the set of eigenvectors K. In 
the present application we select the lowest states obtained in each channel. It turns out 
that these states also have the largest overlap with the model space. Their number is given 
by the number of basis states satisfying 2n + / < iV max . 

Finally, the two-body effective interaction is determined from the two-nucleon effective 
Hamiltonian, obtained from (|9|), as V 2e g = H 2e g — Hq2- Apart from being a function of the 
nucleon number A, V 2e s depends on the HO frequency Q and on the model-space defining 
parameter iV max . It has the important property that V 2c s V12 for iV max — > 00 following 
from the fact that u —>■ for P — > 1. 



E. Three-body effective interaction calculation 

In the standard shell-model calculations the effective interaction is limited to a two-body 
effective interaction. In our approach the limitation to a two-body effective interaction 
means that we use V^g, computed as discussed in the previous subsection, and neglect all 
the three- and higher-body clusters appearing in the expansion (|HJ). The formalism laid 
out in the present paper allows us, however, to employ three-body or even higher-body 
interactions in a straightforward manner. Therefore, we are in a position to go beyond the 
two-body effective interaction approximation. One way to procceed is to evaluate the three- 
body term appearing in Eqs. ([H||l4|,[n|). That term depends, through 5*123 m E<1- (|15|) , on 



the operator u, appearing in Eq. (|17|), used for the construction of the two-body effective 
interaction. It does not, however, guarantee the model space and the Q-space decoupling 
for the three-nucleon system in a similar way as in Eq. ([!(]) for the two-nucleon system. 
To improve on this point, we, therefore, choose a different way of deriving the three-body 
effective interaction. In our previous paper |19| we introduced a procedure applicable for 
nuclei with A > 4 that calculates the three-body correlations directly and guarantees the 
model space and Q-space decoupling for the three-nucleon system. 

In order to calculate the three-body effective interaction in this way, we first re-write the 
two-body interaction term of the Hamiltonian (P||TTD using 

S V *i = S + V* + V jk ) , (19) 

i<j=l i<j<k=l 

valid for A > 3. Then, formally, our Hamiltonian consists only of one-body and three- 
body terms. We now calculate the three-body effective interaction that corresponds to 
Vij + Vik + Vjk from the three-nucleon system condition 

Qse-^ih + h 2 + h 3 + V 12 + V 13 + V 23 )e s(3) P 3 = , (20) 

in complete analogy to Eq. (TTH) . Note that is different from 5*123 an d is determined 
by (pop. The three-body effective interaction is then obtained utilizing the solutions of the 
three-nucleon system solutions from Eqs. (§,[|[|) for the Hamiltonian 

H° = h 1 + h 2 + h 3 + V 12 + Vi 3 + V 23 • (21) 
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As the interaction depends only on the relative positions of nucleons 1, 2 and 3, the three- 
nucleon center of mass can be separated, when solving Schrodinger equation with H%. Sim- 



ilarly as for the two-nucleon Hamiltonian ([TBI), the center-of-mass term is not considered in 



the effective-interaction calculation. We obtain the three-nucleon solutions corresponding 
to the Hamiltonian fl2~l~|) by, first, introducing Jacobi coordinates, as described later in this 
paper and also as described in our previous papers [|H||19|,[2(| and, by, second, the introduc- 



ing for the interactions V12, V13, V23 the two-body effective interactions corresponding to a 
large space characterized by N 3max « 30 and derived according to the procedure described 
in the previous subsection. A space of this size is sufficient for obtaining exact or almost 



exact solutions of the three-nucleon problem [lq|. We note that 

V ij = V N (r i -r j )-—(r i -r j ) 2 , i,j = 1,2,3, (22) 

yields a three-nucleon bound system for A > 4. Therefore, as in the case of the two- 
body effective interaction calculation, the Lee-Suzuki approach is applicable, as described 
in Subsection [11 B|, in a straightforward way. We use the solutions of the three-nucleon 



system to construct the operator to and, then, the three-body effective interaction directly 
from these solutions by applications of the relations (j|) and @. The eigensystem E k , \ k) is 
obtained from the solution of Schrodinger equation (JJ) for the Hamiltonian (pi]) with the 
interactions replaced by V^s,^-, where the two-body effective interaction corresponds to 
the space defined by iV 3max , as discussed above. The model space, defined by the projector 
P 3 , is characterized by the maximal allowed number of HO excitations A" max , where A" max < 
-/V3max- It is spanned by all states satisfying the condition jV 3 < A" max with A" 3 given as 
A" 3 = 2n + I + 2Af + C when basis ( ^9|) is used, or, using the basis (j32|) A" 3 = N . The 
Q-space defined by Q3 is spanned by states with total number of HO excitations iV max < 
^3 < N 3max . Typically, we solve the three-nucleon system and construct the three-body 
effective interaction only for three-nucleon channels with J 3 = l/2 ± ,3/2 ± and T 3 = 1/2. 
For all other channels the two-body effective interaction Yh<j=i ^2eff,ij corresponding to the 
model space characterized by A" max is used instead. 

In order to construct the operator to @ we need to select the set of eigenvectors K. As 
in the two-body effective interaction calculation, the lowest states obtained in each channel 
are selected. Again, those states also have the largest overlap with the model space. Their 
number is given by the number of basis states with the total number of HO excitations A3 < 
Am ax . The three-body effective interaction is determined from the three-nucleon effective 
Hamiltonian, obtained from (||), as V^af = -f^eff — -#03 with H 03 equal to hi + h 2 + h 3 minus 
the three-nucleon center-of-mass HO term. Like the two-body effective interaction V^efr, 
apart from being a function of the nucleon number A, depends on the HO frequency Q 
and on the model-space defining-parameter A" max . In addition, it also depends on the choice 
of A3 max . Obviously, A" 3max must be sufficiently large, in order to make this dependence 
negligible. The limiting properties of are as follows: V^g — > X^<j=i ^2eff,ij for A" max — > 
iV"3max and V 3efi -> E?<j=i Vij for A" max , N 3max -> 00. 



S 



III. TRANSLATION ALLY-INVARIANT HARMONIC-OSCILLATOR BASIS 



As discussed in the previous sections, by using the effective interaction theory we arrive 
at a Hamiltonian that has the following structure 



H AcS 



E 



2m 




(23) 



eff 



with the interaction term depending on relative coordinates (and/or relative momenta) only. 
The center-of-mass dependence appears only in the one-body HO term. Consequently, by 
performing a transformation to Jacobi coordinates, the center-of-mass degrees of freedom 
can be removed. 



A. Jacobi coordinates 



We work in the isospin formalism and consider nucleons with the mass m. A generaliza- 
tion to the proton-neutron formalism with unequal masses for the proton and the neutron is 
straightforward. We will use Jacobi coordinates that are introduced as an orthogonal trans- 
formation of the single-nucleon coordinates. In general, Jacobi coordinates are proportional 
to differences of centers of mass of nucleon subclusters. 

For our purposes we need three different sets of Jacobi coordinates. The first set, i.e., 



— * 

£a-2 



n + r 2 + . . . + r A ] , 




fi + r 2 + . . . + f A _ 2 ) - r A -i 
(fi + r 2 + . . . + r^-i) 



(24a) 
(24b) 
(24c) 

(24d) 
(24e) 



is used for the construction of the antisymmetrized HO basis. Here, £o is proportional to 
the center of mass of the A-nucleon system. On the other hand, £ p is proportional to the 
relative position of the p + 1-st nucleon and the center of mass of the p nucleons. 

Another set that is suitable for the basis expansion, when two-body interaction matrix 

— * — * — * 

elements need to be calculated, is obtained by keeping £o,£i ? • • -Ci-3 an d introducing two 
different variables, as follows 

(25a) 

■ 1 



£o ; £i, • • • £4-3 , 



VA-2 = 



Va-i 




n + r 2 + . . . + r A _ 2 ) 



(25b) 
(25c) 
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Eventually, a set suitable for the basis expansion, when three-body interaction matrix 
elements need to be calculated, is obtained by keeping £ > £1, • • • £a-4 an d Va-i from the 
previous set and introducing other two different variables 



^A-3 = 



3(A-3) 



A-2 



A-3 



(ri + r 2 + 



+ ^4-3) - - (fA-2 + ^4-1 + rx) 



A-2 



r^A-i 



(r*A-i + f A ) - r*A- 



2 ^ 



7"aI 



(26a) 
(26b) 

(26c) 

(26d) 



Identical transformations, as in Eqs. 



are also introduced for the momenta. 



Let us note that the one-body HO potential in Eq. (|23"D transforms as 



2 -2 
r 



™^ 2 5 + E >^ 



A-l 



2^2 
P ' 



(27) 



and the kinetic term transforms in an analogous way. As the interaction does not depend on 
£0 or the center-of-mass momentum, the center-of-mass HO term can be omitted. Moreover, 
we can use the HO basis, depending on coordinates £ p , p — 1, 2 . . . A — 1, e.g., 



A-l 

n Civ?) . 
P =i 



(28) 



for our calculations. First, let us remark that due to the othogonality of the transformations 
the same HO parameter is used for all HO wave functions. Second, we may use 

), as similar relations like 



any of the sets (|2^j25|j26]) for the basis construction of the type ([ 
([Z7|) hold for all the sets. Third, the basis ( [ZBD is not antisymmetrized with respect to the 
exchanges of all nucleon pairs. The antisymmetrization procedure will be discussed in the 
following subsections. 



B. Antisymmetrization for the three- nucleon system 

We discussed the antisymmetrization of the translationally-invariant HO basis for the 
three-nucleon system in our previous papers |l8| , p!9| , p6| . One starts by introducing a basis 
following from Eq. ( p8|) that depends on Jacobi coordinates £1 and £2, defined in Eqs. (p|). 
e.g., 

\{nlsjt;NCJ)JT) . (29) 

Here n, I and J\f, C are the HO quantum numbers corresponding to the harmonic oscilla- 
tors associated with the coordinates (and the corresponding momenta) £1 and £ 2 > respec- 
tively. The quantum numbers s,t,j describe the spin, isospin and angular momentum of 
the relative-coordinate two- nucleon channel of nucleons 1 and 2, while J is the angular 
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momentum of the third nucleoli relative to the center of mass of nucleons 1 and 2. The 
J and T are the total angular momentum and the total isospin, respectively. Note that 
the basis (|29| ) is antisymmetrized with respect to the exchanges of nucleons 1 and 2, as the 
two-nucleon channel quantum numbers are restricted by the condition (— \y+ s+t = — l. It is 
not, however, antisymmetrized with respect to the exchanges of nucleons 1 •*->• 3 and 2^3. 
In order to construct a completely antisymmetrized basis, one needs to obtain eigenvectors 
of the antisymmetrizer 

*=i(l+T«+T«), (30) 

where and T^> are the cyclic and the anticyclic permutation operators, respectively. 
The antisymmetrizer X is a projector satisfying XX = X. When diagonalized in the basis 
(H|), its eigenvectors span two eigenspaces. One, corresponding to the eigenvalue 1, is 
formed by physical, completely antisymmetrized states and the other, corresponding to the 
eigenvalue 0, is formed by spurious states. There are about twice as many spurious states 
as the physical ones PB . 



Due to the antisymmetry with respect to the exchanges 1^2, the matrix elements in 
the basis ( p9|) of the antisymmetrizer X can be evaluated simply as (X) = |(1 — 2(7~ 23 )), 
where T 23 is the permutation corresponding to the exchange of nucleons 2 and 3. Its matrix 
element can be evaluated in a straightforward way, e.g., 



((n 1 hs 1 ] 1 t 1 ;N 1 C 1 J 1 )JT\T 23 \(n 2 l 2 s 2 j 2 t 2 ;Af 2 C 2 J 2 )JT} 




\ h si jx 

ls { L S J 

x{n 1 l 1 M 1 C 1 L\M 2 C 2 n2l2L) 3 , (31) 

where Ni = 2rij + /, + 2J\fi + £i,i = 1, 2; j = + 1; and (riiliNiCiL\N 2 C 2 n 2 l 2 L)2, is the 
general HO bracket for two particles with mass ratio 3 as defined, e.g., in Ref. ||27|| . The 



expression (31) can be derived by examining the action of T 23 on the basis states (|29|) . That 
operator changes the state \nl(^i), J\f£(t; 2 ), L) to |nZ(^' / 1 ) , Af£(!;' 2 ), L), where = 1,2 are 
defined as £j, % = 1,2 but with the single-nucleon indexes 2 and 3 exchanged. The primed 
Jacobi coordinates can be expressed as an orthogonal transformation of the unprimed ones. 
Consequently, the HO wave functions depending on the primed Jacobi coordinates can be 
expressed as an orthogonal transformation of the original HO wave functions. Elements of 
the transformation are the Talmi-Moshinsky HO brackets for two particles with the mass 
ratio d, with d determined from the orthogonal transformation of the coordinates, see, e.g., 



Ref. [27 



The resulting antisymmetrized states can be classified and expanded in terms of the 
original basis (|^) as follows 

\NiJT) = J2(nlsjt;AfCJ\\NiJT}\(nlsjt;AfCJ)JT) , (32) 

where N = 2n + / + 2J\f + L and where we introduced an additional quantum number i that 
distinguishes states with the same set of quantum numbers N, J, T, e.g., i — 1, 2, . . . r with 
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r the total number of antisymmetrized states for a given N, J, T. It can be obtained from 
computing the trace of the antisymmetrizer X [21 



TrX 



NJT 



(33) 



C. Antisymmetrization for the A- nucleon system 

The construction of a translationally-invariant antisymmetrized HO basis for four nucle- 
ons that contains an antisymmetrized three-nucleon subcluster was described in our earlier 
papers fl9y26|j . In this subsection we generalize and simplify this construction to the case of 
an arbitrary number of nucleons A. The starting point is a basis that contains an antisym- 
metrized subcluster of A — 1 nucleons, e.g., 

|(A/ A _ 1 2 A _ 1 J A _ 1 T A -i; tia-iIa-iJa-i)JT) , (34) 

with the (^4 — l)-nucleon antisymmetrized state |iV A _iz A _i J a _iTa-i), depending on the Ja- 
cobi coordinates ^1,^2, • • • £4-2 (p4]) , and the state [ua-iIa-iJ^a-i) that represents the last 
nucleon depending on the Jacobi coordinate For a four nucleon system the state 

|^Va-i«a-i>7a-i^a-i) is identical to the state introduced in Eq. fl32"l) . The basis ([34]) is 
not antisymmetrized with respect to exchanges of the last nucleon with the others. How- 
ever, due to the antisymmetry of the A — 1 nucleon subcluster, the matrix elements of the 
antisymmetrizer X in the basis (S3]) simplifies dramatically, e.g., 



(X) = j[l-(A-1)(T AA ^)] (35) 

with Ta.a-i the transposition operator of the v4-th and the A — 1-st nucleon. Its matrix 
element can be computed in a straightforward way as 

((^a-il^a-il^a-il^a-il; h. a _i]J A -ilJa-il) JT[ 
x^a,a-iI(-^a-ir^a-ir^a-ir2a-ir; ^a-ir^a-irJa-ir)<^2/ 

= 8n l ,N r ^ (-/Va-2^ A-2 Ja~2?A~2 i WA-2]Ja-2Li7a-2l| |^V A -1L«A-1L'7a-1lT A _il) 
X (■^V A -2*A-2^A-27a_2; ^A-2R^A-2Rv7a-2r| I-^A-IR^A-IR^A-IrTa-ir) 

yf, T ( n T A-lL+T A -lR+JA-2L+JA-2R / \ ^A-2 ^A-Ir] 

{ 2 1 J A-1L J 

{^A-2 Ja-2L </a-1L 
i7a-2R K Ja-vl 
Ja-xk Ja-xr J 

Ia-2L Ia-1R K 1 I / A -2R ^A-1L ^11 ^A-1L ^A-2R ^ 1 
> i7a-1R Ja-2X \ J \ Ja-XL Ja-2R 2/I ^A-1R ^A-2L £ J 

xL 2 (-l) iA - 2R+iA - 1L+L (n A _ 1L / A _ 1L n A _ 2L /A-2L ^|n A - 2 R/A-2R^A-lR/A-lR £>A(A-2) , (36) 

where iV x = iV A _i X + 2n A _ iX + Za-ix, 

X = L,R and (n A _i L /A-iL^A-2L^A-2L^|^A-2R^A-2R^A-iR^A-iR^)A(A-2) is the general HO 



x 
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bracket for two particles with mass ratio equal to A(A — 2). The expression (|36|) reduces to 
Eq. (19) in Ref. |l!| for A = 4. Let us stress the important property of the antisymmetrizer 



matrix, namely its diagonality in N A = -ZV L = Ar. Consequently, we may impose a model 
space restriction of the type N A < A" max and still obtain all the antisymmetrized states 
within that model space. 

The antisymmetrized states are obtained by diagonalizing the antisymmetrizer X (|35D 
in the basis fl3"l| ) or, more efficiently, by employing the method introduced in Ref. pSJ that 
does not require us to compute all the matrix elements of the antisymmetrizer. In order to 
get the basis for the A-nucleon system we need to set up an iterative procedure that starts 



with the calculation of the three-nucleon basis (|32| ) and procceeds to four nucleons and so 
on up to A. The resulting antisymmetrized states can be classified and expanded in terms 
of the original basis fl34|) similarly as in the three-nucleon case 



\N A i A JT) = ^(^A-lU-l^A-l^A-l^A-l^A-lJA-lHiVAU^) 

\(N A -iiA-iJA-iTA-i;nA-ilA-iJA-i)JT) , (37) 

where A^a = A^a-i + 2nA-i + Ia-i- In Eq- (0)' ^ ne additional quantum number i A distin- 
guishes states with the same set of quantum numbers N A , J, T, e.g., %a = 1, 2, ... r with r 
the total number of antisymmetrized states for given N A , J,T. Again, it can be obtained 
from computing the trace of the antisymmetrizer X 

r = TtX NaJT . (38) 



D. Recoupling to basis containing two- and three-body clusters 

The A-nucleon antisymmetrized basis, obtained as described in the previous subsection, 
is given as an expansion of the basis (|34]), as shown in Eq. (|37D. In this form, it is not, 
however, in general suitable for calculations with two-body or three-body interactions. In 
order to facilitate the calculations with two-body or three-body interactions we need to 
expand the antisymmetrized states \N A i A JT) in a HO basis consisting of antisymmetrized 
subclusters of A — 2 and 2 nucleons depending on the Jacobi coordinates (^) or consisting 
of antisymmetrized subclusters of A — 3 and 3 nucleons depending on the Jacobi coordinates 
(PB|), respectively. 

For a calculation for an A > 3 system with a two-body interaction, we need the following 
expansion matrix element 

((N A - 2 iA-2JA-2T A -2;(nlsjt,N£)J)JT\N A i A JT) 
= X^a-iU-i^a-iTa-i; ^a-i^a-iJa-iII^aU^) 

X (^A~2^A-2^A-2^A-2; ™A-2^A-2i7a-2 | I^A-l^A-l^A-l^A-l) 
X J A -lJA-2^iA-ljs(-l) JA - 2+JA - 1+JA - 2+J+i+£+S+ ' A - 1+ZA - 2 



XI 



-l)T A -2+T+if A _^ 1 T Af2 | T A1 1 I ^ I ^-a i , j j j 



2 T f j } L I J \ \ ^A-l J J 




XL (nlMC £|«A-l^A-l^A-2^A-2 L) A , (39) 
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where we used an orthogonal transformation of the Jacobi coordinates £,a-2i£,a-i and 
f]A-2,VA-i- In the state | (-^-2^.-2^-2^-2; (nlsjt,J\f£)J^)JT), the antisymmetrized sub- 
cluster |A^a-2*a-2<7a-2^a-2) depends on the Jacobi coordinates £1, £2, • • • £4-3 • The two- 
nucleon channel state \nlsjt) depends on the Jacobi coordinate 7X4—1 and the HO state 
\JVC) that describes the relative motion of the two subclusters is associated with the Jacobi 
coordinate f)A-2, given in Eq. (|25|). 

When this expansion of \N A i A JT) is used a matrix element of a two-body interaction in 
the basis \N A i A JT) can be evaluated in a simple manner, e.g., 

( E v a) = \ A ( A ~ 1)<V(v^a-i)> , (40) 
i<j=l z 

and the matrix element on the right-hand side is diagonal in all quantum numbers of the 
state I (A^a-2^a-2^a-2^a-2; ( n ^sjt, J\f C)J)JT) except n, I for an isospin invariant interaction. 

Similarly, for a calculation for an A > 5 system with a three-body interaction, we need 
the following expansion matrix element 

((N A ^i A ^J A ^T A ^(N 3 i 3 J 3 T 3 ,AfC)J)JT\N A i A JT) 
= E( Ar A-2U-2JA-2T A -2; (nlsjt, N'C')J'\\N A i A JT) 
x(N A 

-3*A-3"^A-3^A-3j ^A-3^A-3i7a-3 | |^ v A-2^A-2^A-2^A-2/ 

x(nlsjt; n A _ 3 l' A _ 3 J A _ 3 \\N 3 i 3 J 3 T 3 } 
xJ A . 3 J' A . 3 JJ'J A _ 2 J 3 (_l)^'+JA-a+J+i-HA-3+ \ 

\- r A-3+T 3 +Trfn rp j ^A-3 2 ±A ~ 2 I J ^ A -3 J 3 

t T T 3 (1 J' Ja-3 Ja-2 

Jk-3 h 3 



xL 2 (-l) L {n' A _ 3 l' A _ 3 NC L\n A - 3 l A -^TC' L) 2(A - 3) . (41) 



Here we used an orthogonal transformation of the Jacobi coordinates £a-3, Va-2 and 
Va-sJa-2 given in Eqs. In the state |(A^a-3U-3 Ja-sTas] (N 3 i 3 J 3 T 3} N C)J)JT) , 

the antisymmetrized subcluster |A^a-3«a-3<7a-37a-3) depends on the Jacobi coordinates 
£1, £2, . . . £a-4- The three-nucleon antisymmetrized subcluster \N 3 i 3 J 3 T 3 ) depends on the 
Jacobi coordinates $a-2,Va-i and the HO state \J\f£) that describes the relative motion of 
the two subclusters is associated with the Jacobi coordinate $a-3, given in Eq. (|26|). 

When this expansion of \N A i A JT) is used a matrix element of a three-body interaction 
in the basis \N A i A JT) can be evaluated in a straightforward way, e.g., 

( J2 V ijh ) = ±A(A-l)(A-2)(V($ A - i ,fr A - 1 )), (42) 
i<j<k=i " 

and the matrix element on the right-hand side is diagonal in all quantum numbers of the 
state |(A^A-3iA-3^A-3^A-3i (N 3 i 3 J 3 T 3 ,J\f C)J)JT) , except N 3 and i 3 , for an isospin invariant 
interaction. 
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E. One-body densities 



In the formalism of the translationally-invariant shell model, one can introduce one- or 
higher-body densities, as in the standard shell-model formulation. Let us consider a general 
one-body operator, e.g., 

d^ = j^d^\f l -R,a l ,T l ). (43) 
i=i 

Its matrix element between antisymmetrized states depending on the Jacobi coordinates 
(Ell) can be written schematically as 



(6™) = A{d^\-^^ 7 ±i A _ ll a A ,T A) ) . (44) 

In a more detailed form, we can express a reduced matrix element between two eigenstates 
of a Hamiltonain corresponding to A-nucleon system as 

(A;EJ T T\\\d {kT) \\\A; E'J'^'T') 

= J2(A;ErT\(N A „ 1 i A _ 1 J A _ 1 T A _ l] n A _ 1 l A _ 1 J A _ 1 )JT) 
x ((iV A -iU-i Ja-iTa-i; n' A _ i r A -iJLi)J'T'\A; ET'T) 

X Jjl(-l)J^+K+J+Jk-, | ^A-l ^ | ^/ ( _ 1) T A _ 1+ r+T + \ | Ta_i | T' | 

x (n A _ 1 Z A _ 1 J A _ 1 1 1 1 AO^(- J ^-^O-i, *a, ta) I IK-A-i^a-i) , (45) 



where we used expansions of the eigenstates in the basis (|37|). Apparently, the matrix 
element (fi5|) factorizes into products of one-body reduced matrix elements and one-body 
densities. 

One can introduce two- or higher-body densities in a similar way. 



IV. RESULTS FOR FEW-NUCLEON SYSTEMS 

We have written a computer code for calculations using the formalism presented in 
Sections [TI| and [Til]. We performed test calculations for few-nucleon systems up to A = 8. 
The code reproduces results obtained using the many-fermion dynamics (MFD) shell-model 
code [0 > when a two-body effective interaction is employed and when the model space does 



not require more than 9 major HO shells, i.e., the limits of the version of the MFD code 
we have. Due to the computational complexity of the antisymmetrization procedure we 
expect that the present formalism can, at the present time and in the current formulation, 
improve significantly on the results obtainable by the MFD code for A < 6 and to some 
extent for A = 7 and A = 8. Let us further remark that although the antisymmetrization is 
complicated it needs to be done only once for given A, N A , J and T. 

In this paper we present results for A = 3 and A = 4 systems, while calculations for 
larger A will be published separately at a later stage. We investigated the A = 3 and A = 4 
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systems in the framework of the present formalism in two previous papers [|18j and |19 
respectively. However, as the newly developed code is more efficient, we were able to extend 
the calculations to larger model spaces. In addition, in this paper we present results obtained 
using the momentum-space dependent non- local CD-Bonn NN potential There are no 



published results for the four-nucleon system interacting by CD-Bonn potential up to now. 

We work in the isospin formalism. As the CD-Bonn NN potential breaks the isospin 
and charge symmetry we construct an isospin invariant potential for the T = 1 two-nucleon 
channels by taking combinations of V np , V pp and V nn . For 3 H we take \V np + |Km, for 3 He 
we take \V np + §Vp P , while for 4 He we use \V np + \V pp + \V nn . The Coulomb potential is 
added to the CD-Bonn V pp potential. Similarly, for the nucleon mass we use the proton and 
neutron mass combinations, e.g., m = -j(Zm p + Nm n ) with Z and N the number of protons 
and neutrons, respectively. 



A. 3 H and 3 He 

For the A = 3 system we use the two-body effective interaction calculated as described in 



Subsection |II D| . As the effective interaction depends on the model-space size, characterized 
by A max , and the HO frequency Q, we investigate the dependence of observables on those 
two parameters. We performed calculations in the model spaces with A max up to 34 for 
a wide range of HO frequencies Q, e.g., hQ = 19 ... 32 MeV. Our ground-state results are 
presented in Figs. |1]-[| In Fig. ^ we show the 3 H ground-state energy dependence on the 
model-space size in the range of iV max = 4 — 34. Different full lines connect results obtained 
with different HO frequencies. The dotted line represents the 34-channel Faddeev equation 
result -8.00 MeV It is apparent that our results converge to the Faddeev equation result 
as A max increases. We note that the fundamental approximation used in our approach is 
the negligence of the three-body clusters in the expansion (|T3|). Such clusters can give both 
positive and negative contribution to the ground-state energy. Our calculation is not a 
variational calculation. Therefore, we cannot expect a convergence from above. As seen 
from Fig. [I] our results converge both from above or below, with some oscillations possible, 
depending on the HO frequency employed. In Fig. [| we present the same as in Fig. [1] for 
A max in the range from 14 to 34 using an expanded energy scale. We can see that, even if 
a complete convergence has not been achieved for A max = 34 in the whole range of Q used, 
it is possible to interpolate from the different curves to obtain the converged result. For 
example, the line corresponding to Ktt = 28 MeV remains almost constant for A max > 18 at 
the value about -8.00 MeV. We present the interpolated results in Table |. Our 3 H result 
-8.002(4) MeV is in a good agreement with the 34-channel Faddeev calculation. We note, 
however, that in our calculations we used all the two-nucleon channels up to j = 6. We 
should, therefore, compare with the result -8.014 MeV obtained by Nogga et al. [|J, where 
all channels with j < 6 were used. Consequently, it appears that we are missing about 10 
keV in binding, most likely due to imprecisions in the two-nucleon systems solutions used 
for construction of the effective interaction. 

In Fig. |3|, we present the same dependence for 3 He as in Fig. [| for 3 H. Our interpolated 
ground-state energy result together with the error estimate is given in Table |. It should 
be noted that the CD-Bonn potential gives a realistic prediction for the binding energy 
difference of 3 H and 3 He, which is experimentaly 0.764 MeV. On the other hand, the absolute 
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value of the binding energy is underestimated compared to the experimental values, 8.482 
MeV for 3 H and 7.718 MeV for 3 He, by about 400 keV. It is, however, only a half of what 
one gets with the local coordinate-space potentials like Nijmegen, Reid or Argonne. 

The model-space size and Q dependence of the 3 H point-nucleon matter rms radius is 
presented in Fig. |j. We observe a convergence and a saturation with the model-space size 
increase. In Table | we show point-proton and point-neutron rms radii for both 3 H and 
3 He extrapolated from the calculations in the largest model spaces together with the error 
estimates. In addition, the calculated magnetic moments are also presented. These can be 
compared to the experimental values of +2.979 for 3 H and —2.128 for 3 He. The 
calculated values were obtained using a one-body Ml operator with bare nucleon g-factors. 



B. 4 He 

The calculations for 4 He were performed in the model spaces up to iV max = 16 in a wide 



range of HO frequencies fi. This is an extension of our previous A = 4 calculations ||19|| , where 
model spaces only up to iV max = 14 were utilized and a narrower range of Q was investigated. 
We performed separate calculations both with two-body effective interactions, computed 
as described in Subsection [11 L)[ and with three-body effective interactions, computed as 
discussed in Subsection [11 E] . The three-body effective interactions were calculated for the 
three-nucleon channels with T 3 = \ \ and J 3 T 3 = | \ using iV 3max = 32 and iV 3max = 28, 
respectively. 

Our ground-state energy results are presented in Figs. § and |j. We investigate the 
dependence on both iV max and Q. The two figures show mostly the same points, plotted in 
the first function of iV max and in the second case as a function of hil. The dotted 

lines connect the results obtained using the two-body effective interactions, while the full 
lines connect the results obtained with the three-body effective interactions. We observe a 
decrease of dependence on both iV max and TiVL as the model-space size increases. In particular, 
the curves in Fig. || become more flat with increasing iV max . Similarly as in our previous 
study ji~9|l , it is apparent that calculations done with the three-body effective interaction 



show weaker dependence on both Vt and iV max and demonstrate faster convergence. Due 
to the higher complexity of those calculations, we present results only for hQ = 13, 22, 31 
and 40 MeV. As can be seen from Figs. [5] and [| for HQ less than about 40 MeV the 
binding energy decreases with increasing model-space size while for larger HQ it begins to 
increase with increasing iV max . Similarly as for the A = 3 system, we are in a position to 
interpolate the converged ground-state energy result, though with a lower accuracy. Based 
on the results presented in Figs. [5] and || we estimate the CD-Bonn 4 He ground-state energy 
to be -26.4(2) MeV. The experimental binding energy of 4 He is 28.296 MeV. The CD-Bonn 
thus underbinds 4 He by about 2 MeV. It is again only about a half of underbinding that 
one gets with, e.g., Argonne V18 with the calculated 4 He binding energy 24.1 MeV [pi]]. 



In our approach we obtain the ground state as well as the excited states by diagonal- 
izing the Hamiltonian. In Figs. [7] and |8|, we present the model-space-size dependence of 
both the ground-state and the first excited + state energies and point-nucleon rms radii, 
respectively, obtained in calculations with the three-body effective interactions. Compared 
to the ground state, we observe a much stronger dependence by the excited-state energy 
and nucleon rms radius on both Q and -/V max . The significantly different convergence rate 
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of the ground state and of the first excited + state manifests the different nature of the 
two states. A possible interpretation of this observation is that the excited + state is 
associated with a radial excitation and, thus, it is more sensitive to the HO basis used in 
our calculations. Although we cannot extrapolate the energy or point-nucleon radius of this 
state, our calculations show that its excitation energy is below 21 MeV and the radius is 
larger than 3 fm. 

Our ground-state results are summarized in Table |. 



C. Charge form factors 



A sensitive test of the wave-functions obtained in our calculations is the evaluation of 
form factors. In this subsection we compare the charge form factors obtained with the CD- 
Bonn wave functions and those obtained in an identical calculation with the Argonne V8' 
NN potential defined in Ref. ]10j- We note that a similar comparison was performed by Kim 
et al. in Ref. |35| for the Bonn OBEPQ and the Reid NN potentials. In that paper, it was 
found that the 3 H and 3 He form factors differ for the two NN potentials. 

Using the formalism of Ref. ||32|| , we calculated the charge EM form factors and ratio 
of charge strangeness and EM form factors in the impulse approximation. The one-body 
contribution to the charge operator is given by Eq. (15) in Ref. 
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Li. is the fc-th nucleon orbital momentum, 



T 



and 



(46) 



t) are the 



one-body electric and magnetic form factors, respectively. The superscript (a) refers to (p) 
and (n) for proton and neutron EM form factor, respectively, or to (s) for the strangeness 
form factor. We use the parametrization of the one-body form factors as discussed in Ref. 
32| . We note that the one-body strangeness form factors depend on the strangeness radius 
p s , for which we take the value p s = —2.0 as in Ref. |32| and on the strangeness magnetic 
moment p s . Limits on these parameters are to be determined in the experiments at the 
Thomas Jefferson Accelerator Facility (TJNAF). The first strangeness magnetic-moment 
measurement was reported recently pi and an experimental value p s = +0.23, obtained 



with a large error. We use this value in our calculations. 

The elastic EM charge form factors of 3 H and 3 He are presented in Figs. |9] and [IT], 
respectively. We observe a large sensitivity to the choice of the NN potential. Let us 
remark that we used the wave functions obtained in the model space with iV max = 34. We 
investigated the dependence of the form factors on both iV max (= 30,32,34) and HQ and 
found that the dependence is below the resolution of the figures. Our Argonne V8' results 
compare well with those obtained using the Argonne V18 in the impulse approximation 
presented, e.g., in Ref. [01 1. We note that the experimental position of the minima as about 



3.6 fm -1 and 3.2 fm -1 for 3 H and 3 He, respectively, see Ref. and the references therein. 

The 4 He elastic EM charge form factor and the EM charge form factor corresponding to 
the transition to the first excited + state are presented in Figs. [XT] and [H| respectively. 
The three-body effective interactions were used and the iV max = 16 model spaces. For the 
Argonne V8', the current results can be compared to those presented in Ref. [19| obtained 
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using N max = 14. We note that a second minimum appears in our calculated charge form 



factors in a similar position as found in the VMC calculations presented in Ref. ||34j| . The 
elastic charge form factor sensitivity to both iV max and hQ is weak for q below the secondary 
maximum but it increases for larger q. As to the inelastic form factor, there the sensitivity 
is significantly stronger. We believe that the form factors presented in Fig. O are more 



realistic than the inelastic ones given in Ref. [19| that we obtained using iV max = 14. 



In general, for all 3 H, 3 He and 4 He the CD-Bonn results are further from the experimental 
data points than the results obtained using the Argonne V8'. This is in full agreement with 
the calculations in Ref. |£J for 3 H and 3 He. However, in order to make any conclusion about 
the superiority of any of the potentials, one needs to calculate the meson exchange current 
contributions. 

Finally, in Fig. [13| we present the ratio of the 4 He charge strangeness and EM form 
factors calculated in the impulse approximation using both the Argonne V8' and the CD- 
Bonn NN potentials. The ratio of the elastic charge form factors is particularly interesting, 
as it can be experimentaly obtained from the measurement of the parity-violating left-right 
asymmetry for scattering of polarized electrons from a 4 He target. Experiments of this type 
are now under preparation at TJNAF. 



V. CONCLUSIONS 

We presented a translationally invariant formulation of the no-core shell-model approach 
for few-nucleon systems and introduced a general method of antisymmetrization of a HO 
basis depending on Jacobi coordinates. The latter procedure starts with the construction of 
the antisymmetrized basis for three nucleons, then procceeds to four and so on. We derived 
an iterative algebraic formula for computing the antisymmetrized basis for A nucleons from 
the antisymmetrized basis for A — 1 nucleons. In addition, we discussed how to transform 
the antisymmetrized states to bases containing different antisymmetrized subclusters of nu- 
cleons. The chosen approach has the advantage that the antisymmetrizer is very simple and 
that the dimensions of the starting basis, formed by the A — 1 nucleon antisymmetrized sub- 
cluster and the last nucleon, are the lowest compared to bases with different subclustering. 

There are two main advantages of the use of a translationally- invariant basis. First, it 
allows us to employ larger model spaces than in traditional shell-model calculations. Second, 
in addition to two-body effective interactions, three- or higher-body effective interactions 
as well as real three-body interactions can be utilized. The use of higher-body effective 
interactions reduces the dependence on the HO frequency and speeds up the convergence of 
our approach. 

As the antisymmetrization procedure is computationally involved, the practical applica- 
bility of the formalism is limited to light nuclei. In the present formulation we expect that 
significant improvement over the traditional shell-model results can be achieved for A < 6 
and to some extent for A = 7 and A = 8. However, different paths of antisymmetrization 
than that realized here can be chosen for more complex nuclei. On the other hand, there 
exists a sophisticated approach for antisymmetrization of hyperspherical functions depend- 
ing on Jacobi coordinates developed by Barnea and Novoselsky J55|. It makes use of the 
orthogonal group transformations of Jacobi coordinates. That approach can be adapted for 
the HO functions and should lead to a more efficient antisymmetrization. Another issue is 
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the transformation to A — 2 plus 2 and A — 3 plus 3 clusters as discussed in Subsection [IIID 



In principle, it can be avoided by using the reverse ordering of the Jacobi coordinates (P4|) 
|37| . Then, however, the computation of the one-body densities, as described in Subsection 
HIE], would become very difficult. 

We applied the formalism to solve three and four nucleon systems interacting by the 
CD-Bonn NN potential in model spaces that included up to 34hQ and 16M7 HO excitations, 
respectively. For the three-nucleon system our method leads to the exact solution and 
our results are in agreement with the calculations by other methods. For the four-nucleon 
system, we were able to interpolate the ground-state energy solution from the model-space 
and HO-frequency dependencies. Our result with the error estimate is -26.4(2) MeV. There 
have not been any published results by other methods so far for the A = 4 system interacting 
by the CD-Bonn NN potential to which we could compare. In addition to energies, rms radii 
and magnetic moments, we also compared charge form factors obtained using the CD-Bonn 
and Argonne V8' NN potentials. We found a substantial sensitivity to the choice of the 
potential in agreement with results published in Ref. |35| calculated with similar potentials. 

We believe that the method discussed in this paper has the potential to solve the few- 
nucleon problem beyond A = 4. The convergence can still be improved by employing 
four- or higher-body effective interactions in a similar fashion as we used the three-body 
effective interaction. Also the antisymmetrization procedure can be made more efficient. 
Our method has the advantage, compared to, e.g., the GFMC method, that we solve the 
Schrodinger equation by diagonalization. Consequently, wave functions are obtained and 
excited states with identical quantum numbers as the ground state are computed. It would 
be also interesting to compare the present method with a new approach, relying also on the 
HO basis, in which the effective interaction is constructed by solving the Bloch-Horowitz 
equation |3"8fl . Calculations are now under way for A = 5 and A = 6 using the present 
formalism. 
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FIGURES 



FIG. 1. The dependence of the 3 H ground-state energy, in MeV, on the maximal number of 
HO excitations allowed in the model space in the range from N max = 4 to -/V max = 34. The 
two-body effective interaction utilized was derived from the CD-Bonn NN potential. Results for 
HQ = 19, 22, 24, 26, 28, 30 and 32 MeV are presented. The dotted line represents the exact result 
of -8.00 MeV from a 34-channel Faddeev-equation calculation |2l| ]. 



FIG. 2. The dependence of the 3 H ground-state energy, in MeV, on the maximal number of 
HO excitations allowed in the model space in the range from -/V max = 14 to iV max = 34. The 
two-body effective interaction utilized was derived from the CD-Bonn NN potential. Results for 
hfl = 19, 22, 24, 26, 28, 30 and 32 MeV are presented. The dotted line represents the exact result 
of -8.00 MeV from a 34-channel Faddeev-equation calculation |21]. 



FIG. 3. The dependence of the 3 He ground-state energy, in MeV, on the maximal number of 
HO excitations allowed in the model space in the range from -/V max = 14 to ./V max = 34. The 
two-body effective interaction utilized was derived from the CD-Bonn NN potential. Results for 
fin = 19, 22, 24, 26, 28 and 30 MeV are presented. 



FIG. 4. The dependence of the 3 H point-nucleon matter radius, in fm, on the maximal number 
of HO excitations allowed in the model space in the range from iV max = 4 to iV max = 34. The 
two-body effective interaction utilized was derived from the CD-Bonn NN potential. Results for 
htt = 19, 22, 24, 26, 28, 30 and 32 MeV are presented. 



FIG. 5. The dependence of the 4 He ground-state energy, in MeV, on the maximal number 
of HO excitations allowed in the model space in the range from iV max = 8 to iV max = 16. The 
two-body (dotted lines) and three-body (full lines) effective interactions utilized were derived from 
the CD-Bonn NN potential. Results for HQ = 19,22,31,37,40 and 43 MeV are presented. 



FIG. 6. The dependence of the 4 He ground-state energy, in MeV, on the HO energy in the 
range from TiVt = 13 to UQ. = 43. The two-body (dotted lines) and three-body (full lines) effective 
interactions utilized were derived from the CD-Bonn NN potential. Results for -/V max = 8, 10, 12, 14 
and 16 are presented. 



FIG. 7. The dependence of the 4 He ground-state and the first-excited + state energies, in 
MeV, on the maximal number of HO excitations allowed in the model space in the range from 
JV max = 8 to -/V max = 16. The three-body effective interaction utilized was derived from the 
CD-Bonn NN potential. Results for HQ = 13, 22, 31 and 40 MeV are presented. 
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FIG. 8. The dependence of the 4 He ground-state and the first-excited + state point-nucleon 
rms radius, in fm, on the maximal number of HO excitations allowed in the model space in the 
range from iV max = 8 to iV max = 16. The three-body effective interaction utilized was derived from 
the CD-Bonn NN potential. Results for htt = 13, 22, 31 and 40 MeV are presented. 

FIG. 9. The elastic EM charge form factor of 3 H calculated in the impulse approximation. 
Results obtained using the Argonne V8' (dotted line) and CD-Bonn (full line) NN potentials are 
compared. 

FIG. 10. The elastic EM charge form factor of 3 He calculated in the impulse approximation. 
Results obtained using the Argonne V8' (dotted line) and CD-Bonn (full line) NN potentials are 
compared. 

FIG. 11. The elastic EM charge form factor of 4 He calculated in the impulse approximation. 
Results obtained using the Argonne V8' (dotted line) and CD-Bonn (full line) NN potentials are 
compared. The calculations were performed using three-body effective interaction in the model 
space characterized by iV max = 16 and fiQ. = 22 MeV. 

FIG. 12. The EM charge form factor of 4 He corresponding to the transition to the first excited 
+ state calculated in the impulse approximation. Results obtained using the Argonne V8' (dotted 
line) and CD-Bonn (full line) NN potentials are compared. The calculations were performed using 
three-body effective interaction in the model space characterized by N max = 16 and HQ = 22 MeV. 

FIG. 13. The ratio of elastic strangeness and EM charge form factor of 4 He calculated in the 
impulse approximation. Results obtained using the Argonne V8' (dotted line) and CD-Bonn (full 
line) NN potentials are compared. The calculations were performed using three-body effective 
interaction in the model space characterized by iV max = 16 and HQ = 22 MeV. Values of the 
strangeness radius p s = —2.0 and the strangeness magnetic moment [x s = 0.23 were employed. 
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TABLES 



CD-Bonn NN potential 


E gs [MeV] 




V^I) [fm] 


A 4 M 


3 H 


-8.002(4) 


1.608(4) 


1.760(6) 


2.612 


3 He 


-7.248(4) 


1.802(6) 


1.635(4) 


-1.779 


4 He 


-26.4(2) 


1.445(5) 


1.445(5) 





TABLE I. Results for the ground-state energies, point-proton and point-neutron rms radii, and 
magnetic moments obtained for 3 H, 3 He and 4 He using the CD-Bonn NN potential are presented. 
Shown values are based on the results calculated in the largest model spaces used in the present 
study, iV max = 34 for 3 H, 3 He, and N max = 16 for 4 He, respectively. The errors were estimated 
from the dependences on the HO frequency O and on the model-space size characterized by iV max . 
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Q> 


-7.8- 






-8.0- 


UJ 




-8.2- 




-8.4- 




-8.6- 




-8.8- 
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exact 

hn=19 MeV 
hn=22 MeV 
h^=24 MeV 
hn=26 MeV 
hn=28 MeV 
hH=30 MeV 
hn=32 MeV 



At 



max 



Q) -8.0 - 



-8.1- 



-8.3 



1 1 1 1 1 1 1 


i i i i 

3 H -: 


; ci 




)-Bonn : 

i i i i 



12 14 16 18 20 22 24 26 28 30 32 34 36 



N. 



exact 

ha=19 MeV 
ha =22 MeV 
ha =24 MeV 
ha =26 MeV 
ha =28 MeV 
ha =30 MeV 
ha =32 MeV 



max 



■7.1 



■7.2- 



■7.3- 



■7.4- 



t 1 1 1 1 1 1 1 1 1 r 




.751 1 1 1 1 1 1 1 1 1 1 1 

' 12 14 16 18 20 22 24 26 28 30 32 34 36 



hn=19 MeV 
hH=22 MeV 
hH=24 MeV 
hH=26 MeV 
hH=28 MeV 
hH=30 MeV 



N, 



max 



V 



1.7- 



1.6- 



1.5- 



1.4- 



i — i — i — i — i — i — i — i — i — i — i — i — i — i — i — r 




hft=19 MeV 
ha=22 MeV 
ha=24 MeV 
ha=26 MeV 
ha=28 MeV 
hft=30 MeV 
hft=32 MeV 



I ^LJ I I I I I I I I I I I I I I I 

4 6 8 10 12 14 16 18 20 22 24 26 28 30 32 34 36 



N, 



max 




q[fm" 1 ] 




max 



-25 n 1 1 1 1 1 1 1 r 




_31 u i i i i i i i i i i_ 

13 16 19 22 25 28 31 34 37 40 43 



hQ [MeV] 



Nmax = 


16 


Nmax = 


16 


Nmax = 


14 


Nmax = 


12 


Nmax = 


10 


^max = 


8 



V 




max 



